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We report a solution of the tilted-axis cranked HFB equation for ^''''Hf, which shows wobbling 
motion coupled to gamma vibration at high spin (J ~ 60?i). Possible anharmonicity and splitting 
of energy levels are also discussed as a consequence of the wobbling motion with large amplitude. 
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In understanding the excited states of a quantum many-body system vifith finite degrees of freedom, the concept of 

dynamic collective modes has been playing a leading role. For example, in nuclear physics, shape vibrations such as 7 

and /3 vibrations [|l[ , octupole vibrations in super-deformed states ^ , and giant dipole resonances Q are well known. 

In a dilute atomic gas in a Bose-Einstein condensate, which is also a finite system, the dynamic motion displays the 

p^ scissors mode, which is an out-of-phase vibration among states in the super and normal fluidity yj. 

I As well as these vibrations, rotations have been also considered and studied in detail. As a result, rotational bands 

^ ] are observed in nuclei, as in atomic and molecular systems |^,|lj. The relevant rotational forms are essentially uniform 

' ' and one-dimensional, and they are sometimes called "static rotation" because the rotation axis is fixed. 
00 [ Advances in experimental techniques these days allow us to obtain data for excited levels built on very high 

'""^ ' spin states. For understanding these levels, a concept of "dynamic" nuclear rotation is introduced and recently has 
attracted much interest. Indeed, the concept is helpful to comprehend several novel phenomena, such as nuclear 
^ ' wobbling motions in triaxial deformed nuclei pUq], dynamical coupling modes (wobbling modes) between high- and 
T-H ' \ow-K bands M, and chiral vibrations through quantum tunnelling between right- and left-handed chiral rotating 
10 states S. These phenomena have recently been studied from both experimental and theoretical approaches. 

Dynamic rotation is usually considered in terms of the classical mechanics as time-dependent and three-dimensional 

rotation, such as the precession in the rotation of the earth M, or motions of a tippe top in which the direction of 

the rotation axis varies time-dependently in the intrinsic frame pO| . The discussion of these rotations in nuclear 

(^ ■ physics was first presented by Bohr and Mottelson in the context of wobbling motion H] . They considered the motion 

as an analogy of the small orientation fluctuations of a triaxially deformed classical rigid rotor from its static and 

one-dimensional rotation about the principal axis of the inertia tensor. They quantize the Hamiltonian of the rotor 

1— H I under the condition |/i| ~ |/| >> 1, by using the technique similar to the second quantization of the harmonic 

^ ■ oscillator. (/ is the total angular momentum and Ji is the x component of the angular momentum vector.) It is 

(-H then shown that the excited levels corresponding to the nuclear wobbling motion have a vibrational character on 

JV the rotational band, E{I,n^) = ^^t. ' + Tiuj^ [n^ -\- 1/2), where fiuj^ is calculated by the moments of inertia of the 

.^H rotor and J7i is the moment of inertia around the a;-axis. The oscillation quanta n^ take the values n^ = 0, 1,2, • • •. 

J^ In spite of their theoretical prediction, it had been difficult to identify the wobbling states in experiments until the 

JH \ experimental evidence recently reported in ^®^Lu with triaxial super-deformation ||^. With this discovery, there has 

■ - - I been established a new research field in high-spin nuclear physics. 

Nuclear rotation is closely related to the shape of a nucleus or its symmetries [y, so that exotic rotations can 
be seen for exotic shapes. Nuclei possessing axial symmetry follow the static one-dimensional rotation, such as in 
the ground-state rotational band, or g-band. These rotations are low-lying collective modes which are caused by 
the spontaneous symmetry breaking mechanism for the rotational symmetry, that is, the SU(2) symmetry is broken 
into the subgroup U(l) through the nuclear deformation. As the U(l) symmetry is further broken into its discrete 
sub-symmetries, more variations can be expected for the nuclear rotation. For instance, dynamic and three- (or two-) 
dimensional rotations are predicted to emerge in nuclei with triaxial symmetry (a discrete symmetry for a tt rotation 
around any axis, denoted as D2). 

One theoretical method used for the analysis of the wobbling motion is based on the particle-rotor model (PRM) 
p,n2lJl3]. Although the model allows quantum mechanical studies, the assumption of a "rotor" brings a macroscopic 
treatment into the analysis. The microscopic formulations for the dynamic rotations were derived in the framework 
of the random phase approximation (RPA) jlj] , and the time-dependent variational method |lq,f6| . Based on these 
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and all subsequent works ||T^,|l8|,0 , we are now able to study the dynamic nuclear rotations from a microscopic point 



of view. The combination of the tilted-axis cranking model |17 19|-|21j and the generator coordinate method (GCM) 
p3,HM| is one of the methods which are very powerful and useful in numerical sutdies. 

Despite the fact that the experimental discovery of the wobbling mode was made in an odd-mass nucleus (^^'^Lu) 
and the theoretical investigation that the wobbling motion is favoured in odd-A mass nuclei [0| , it is essentially the 
"rotor" (in the PRM) that wobbles. Hence, it is necessary to study the wobbling motion in even-even nuclei ( i.e., 
nuclei consisting of the even numbers of protons and neutrons, where all the nucleons are coupled in the ground state 
through the BCS-type pairing interaction). According to Ref. ||ll[, it is calculated that the yrast state (the lowest 
energy state for a given spin) in 72''Hf becomes triaxial super-deformed (TSD) at / ~ AOTi, having gamma deformation 
7 ~ —20°. (Our convention for the gamma deformation parameter ,7, is defined as in p. 7 of Ref. pj.) Thus, ^^"'Hf is 
expected to show the wobbling mode as a result of its triaxially deformed states at high spin. 

In this paper, we thus investigate the nuclear wobbling motion for ^^^Hf, by means of self-consistent calculations 
based on the tilted-axis cranking model. However, before going forward to the microscopic analysis, it is worth 
reviewing rotation of a classical rigid body because the analogy with classical mechanics is quite helpful to understand 
nuclear wobbling motions. 

Rotations of a classical rigid body can be studied through the corresponding angular momentum vector ||9|,|l2|,E3| . 
If the intrinsic coordinates are chosen so as to diagonalize the inertia tensor, cross sections between a sphere, M — 

Mi+M2+M^ (angular momentum conservation), and an ellipsoid, E — 2j^ + T7" + T7" (^^i^i'Sy conservation), give the 
trajectory of the vector on the sphere. (Mj denotes the i-th component of the "classical" angular momentum.) Fig|^(a) 
shows examples of such trajectories. When the rigid body has axial symmetry (say, along the z-axis), it is well known 
that the trajectories follow "precession" , which means that motions of the vector are circular around the z-axis and 
in the plane parallel to the x — y plane. However, in the case of a triaxial rigid body (Ji ^ J2 ^ Js), the trajectories 
deviate from the precession in two ways: (i) distortion from the circle and (ii) deviation from two-dimensional to 
three-dimensional movement. These deviations imply that additional vibrations are induced by breaking the axial 
symmetry of the rigid body. These vibrational motions in rotations are called "wobbling motions" |2^ . It should 
be noted that these deviations from the precession are small when the angular momentum is polarized near the axes 
{x- or z-axes in Fig.nh . This type of small fluctuation corresponds to the nuclear wobbling motion in Ref. m if the 
total angular momentum |M| is large enough ;|M| ~ \Mi\ » 1. Dashed lines in Fig.|l|(a) correspond to special 
trajectories called "separatrices" , which divide the area into four topologically different domains. The existence of 
such domains is evidence of non-linearity in the system. In fact, in Fig.Wa), there are two separatrices separated at 
points V (Ml, M2, M3) = (0, ±1, 0), and they do not cross each other, as shown in Fig.|l|(b). This subset figure shows 
the time-dependence of the angular momentum vector for the separatrix in the domain z > 0, which is obtained by 
solving the Euler equations (-^ = -^-j—^ujjuik, ) with proper initial conditions. (Indices {i,j, k) should be taken in a 
cyclic manner). Most of the time, the vector stays at V and the motion between the two points is comparatively fast 
once the vector is apart from these points. The other separatrix has similar motion except the sign of M^. 

Now, let us begin the microscopic analysis by means of the tilted-axis cranked HFB method. Our Hamiltonian 
consists of two terms: a spherical part ( the spherical Nilsson Hamiltonian) and a residual part (pairing-plus-Q-Q force). 
The model space and force parameters are chosen in the same manner as in our previous work [ p4[ . The HFB equation 
is solved by means of the method of steepest descent, under constraints on the angular momentum and the particle 
number. Particularly, the angular momentum constraints are expressed as, (Ji) = Jcosip, (J2) = 0, (J3) — Jsinip, 
where f denotes the tilt angle. For details about the method, see Ref. ||l^. An advantage of the present method 
is that the physical quantities like the deformation parameters are calculated self-consistently. It is thus possible to 
describe the deformations as functions of the tilt angle, for instance, 7 = 7(1/?). 

First of all, let us see briefly the results obtained by the principal-axis (one-dimensional) cranking calculations 
{(f — 0°). In the present calculations, there are no super-deformed solutions appearing in the yrast band throughout 
the range < J ^ 807i. This can be explained by that fact that our model space is not large enough to reproduce 
super-deformation. Nevertheless, at J ~ 607i, the yrast states take on gamma deformation (7 ~ —20°,/? ~ 0.13). For 
the purpose of studying the wobbling motion theoretically, emergence of triaxility is enough. The gamma deformed 
states appear after of the gap energy vanishes at J ~ 50h (207i) for protons (neutrons). In the present calculations, 
major contributions to the gamma deformation are alignments of ii3/2 and hii/2 protons in addition to neutron 
alignments at lower spins (J ~ 20h); ii3/2, f7/2 7 and hg/2- 

Next, let us see the results in the tilted-axis (two-dimensional) cranking calculations. The tilted-axis cranked HFB 
states are created from the principal-axis cranked HFB states {ip = 0°) a,t J = 6QT1 as the initial states. The range of 
the gamma deformation is self-consistently determined and given in the range 180° ^ 7 ^ 300°. Due to the symmetry 
in the definition of 7, the range may be converted into —60° < 7 < 60°. The conversion is equivalent to renaming 



the axes in a cyclic manner, for instance, "a;-axis -^ y-axis". As a result, the wobbling motion in this study is in the 
x-y plane. In other words, it can be said that the tilt angle in the present work corresponds to the aziniuthal angle 
(if) while in the previous works [0J2j| it corresponds to the polar angle (9), if the z-axis is chosen for the direction of 
elongation of the nucleus. 

Figs.^(a) and (b) show, respectively, changes of gamma deformation, 7((/?), and energy, E{ip). From (a), couphng 
of the tilting and gamma degrees of freedom is clearly seen. It should be noted that the curvatures of 7(9?) are flat 
around t^p ~ 0° and ±90°. This result implies that the assumption in the PRM studies (fixing the gamma values) pjlq] 
is reasonable as long as the wobbling amplitude is small. For large wobbling amplitude, it seems that the 7 vibration 
is cooperatively induced by the fluctuation of the orientation of the rotation axis, that is, the wobbling motion. This 
situation is schematically depicted in Fig.0(c) . (The /3 values are almost constant, 0.130 ^ /3 ^ 0.135, for the entire 
range of ip). The nuclear shape becomes prolate when the tilt angle reaches ip = ±45°, which correspond to unstable 
extrema in FigJ3(b). At these points, the ip range is divided into several domains from a dynamical point of view. 
Let us call the domains for —45° < (p < 45° domain-I, in which the wobbling motion is around the j/-axis, while the 
other for 45° < (p < 135° domain-II, in which the wobbling is around the x-axis. Domain-I and domain-II are related 
by a discrete symmetry operation called C4Z ( a symmetry for ^-rotation about the z-axis), due to the dynamical 
coupling between ip and 7. The wobbling motions are topologically classified through these domains. Therefore, it 
can be said that the rotations a,t ip — ±45°, ±135° correspond to "separatrices" in classical mechanics. As seen in 
Fig.g(b), the energy curve looks like a harmonic potential near </? ~ 0° and ~ 90°. From the present calculation, it 
is obtained as V{ip) ~ atp^ where a = 1.64 (MeV). As ip approaches the value for the separatrix, the curve deviates 
from the harmonic shape. Similar features are seen in the wobbling motion of a classical rigid body with triaxiality, 
that is, the deviation from planar motion becomes larger as the trajectory reaches a separatrix. 

Let us discuss the possible manifestation of the excited structures of the wobbling motion. In Ref . |]1| , the wobbling 
excitations are given as a harmonic structure like 7ia'w('T-w + 1/2). Although the harmonic levels are expected to 
be seen for the wobbling motion with small wobbling amplitude, i.e., small n^, anharmonicity can be expected for 
wobbling motions with large amplitude, i.e., large Uw, according to the shape of the energy curve. To see the excited 
structure schematically, we solve the one-dimensional Schrodinger equation for a wobbling phonon, i.e., the quantized 
wobbling motion. It is assumed that the tilt angle can be treated as a dynamical variable, and that E{ip) corresponds 
to the potential V{ip). By putting walls with infinite height at ip = ±45°, the wobbling phonon is confined in the 
domain-I for the sake of simplicity. The formula (4-304) in Ref. m for the wobbling phonon energy {Tiuj^) and the 
coefficient of the approximated harmonic potential, a = Mw^^/^j which is obtained earlier, are used to obtain the 
"mass" for the wobbling phonon, M^/fi ~ a/{hijj^)'^ . The expression for the moment of inertia is employed from 

the expression for a rigid body with quadrupole deformation Ji = ^^'^ ( 1 — J -^(5 cos (7 — ^z) ), in which M is 

the mass of the nucleus (1.53 x 10^ MeV [^) and R = 1.2A^/^ (fm) with the mass number A = 164. The calculated 
wobbling phonon energy ?iWw is given to be 70 (keV) where the quadrupole deformations are {(3,j) — (0.13,-20°) 
at iy9 ~ 0° and J = 60?i. The results are shown in Fig. 3. Because V{ip) keeps the harmonic shape fairly well up 
to \ip\ < 0.4 (i.e., 23°), it is seen that the quantized energies for the wobbling motion {n„ < 5) possess the similar 
structure to the harmonic case, i.e., equal energy spacing. On the other hand, for n^ > 5, the structure starts to 
deviate from the harmonic structure. 

Unlike the classical wobbling motion of a triaxial rigid body, tunnelling effects should be taken into account in the 
nuclear wobbling motion. In other words, each wobbling mode classified by domains can jump from one domain to 
another over the separatrices. As a result of the coupling between ip and 7 degrees of freedom, E{ip) has discrete 
6*42 symmetry, implying that a wobbling phonon moves in the four potential wells. It is well known that the double 
well potentials cause the splitting in energy levels due to the parity symmetry in the corresponding Hamiltonian P6[ . 
The states are then expressed by symmetric and antisymmetric linear combinations of the localized states in each 
potential. In the present case, there would be four degenerate states localized in each potential, if the barrier height 
were infinitely high. Let us denote these localized states as |L1), |L2), |R2), and |R1), for the further left potential to 
the further right. For the parity operation tt, "R" and "L" are exchanged. For instance, 7r|Ll) = |R1). In the case of 
finite barrier height, these four localized states are superposed through the tunnelling to produce two symmetric and 
two antisymmetric linear combinations. These states might form a quartet, and be written as, 

|Sl) = i(|Ll) + |L2> + |R2> + |Rl)) (1) 

|S2) = 1(|L1)-|L2>-|R2> + |R1)) (2) 

1A1) = 1(|L1) + |L2>-|R2>-|R1)) (3) 



|A2)=.i(|Ll)-|L2) + |R2)-|Rl)) (4) 

"S" and "A" denote symmetric and antisymmetric states. If each localized state is in the ground state, the nodal 
numbers for the superposed states are 0,1,2, and 3 for "SI", "Al", "S2" and "A2", respectively. The order of the 
corresponding energy levels is thus expected to be Esi < Eai < Es2 < Ea2- 

In summary, we have investigated the possible wobbling motion in ^^"^Hf by means of the tilted-axis cranked HFB 
states. The calculations show that the nucleus has sizeable gamma deformation at high spin (J 2± 60h), and wobbling 
motion coupled to the gamma vibration can be expected. A possible anharmonicity in the wobbling excitations is 
predicted for higher wobbling excitations, based on the calculated energy curve with respect to the tilt angle. In 
addition, a brief discussion is presented to show that a "wobbling quartet" can be expected for the nuclear wobbling 
motion coupled to the gamma vibration at high spin, as a consequence of the tunnelling effect. 
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FIG. 1. (a): Trajectories of the angular momentum vector of a classical rigid body (with a quadrupole shape; 
(/3,7) = (0.4,-20°)). Separatrices are denoted by chains, (b): Time dependence of angular momentum components on 
the separatrix in the region z > (written as "Separatrix A (2 > 0)" in the figure). 
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FIG. 2. (a) Gamma deformation with respect to tilt angle, 7(1/'); (b) Energy curve with respect to tilt angle, E{ip); and (c) 
a schematic picture of wobbling motion coupled with gamma vibration. At (/3 = ±45° the nucleus has the axial symmetry, and 
the state corresponds to the separatrix in classical mechanics. 
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FIG. 3. Excited levels of the wobbling phonon. Left: the solid curve is the energy (corresponding to Domain-I) calculated by 
the tilted-ajds cranked HFB method for J = QQfi, while the dashed line is the harmonic approximation of the above line. The 
range for the tilt angle corresponds to |</?| < 45°. Horizontal lines denote the quantized levels for the wobbling motion. Right: 
A simple line indicates fiUL!„{n^ + 1/2), while the asterisks connected by lines show the quantized wobbling energies calculated 
in this study. 



